Driven dynamics and rotary echo of a qubit tunably coupled to a harmonic oscillator 
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We have investigated the driven dynamics of a superconducting flux qubit that is tunably coupled 
to a microwave resonator. We find that the qubit experiences an oscillating field mediated by off- 
resonant driving of the resonator, leading to strong modifications of the qubit Rabi frequency. This 
opens an additional noise channel, and we find that low-frequency noise in the coupling parameter 
causes a reduction of the coherence time during driven evolution. The noise can be mitigated with 
the rotary-echo pulse sequence, which, for driven systems, is analogous to the Hahn-echo sequence. 



Circuit quantum electrodynamics implemented with 
superconducting artificial atoms and microwave res- 
onators has emerged as a framework for studying on-chip 
light-matter interactions [1-3 . It has enabled a range 
of experiments including lasing [3], the creation [5HZ] 
and detection [8] of arbitrary Fock states, and microwave 
photon-correlation measurements [9 , 10 . Microwave res- 
onators also provide a means to couple distant qubits 
[TTJ [12] and, in this role, have been used to implement 
quantum algorithms in superconducting circuits [13] and 
to develop quantum computer architectures [14]. How- 
ever, the coupling of a qubit to a resonator also influences 
the qubit coherence, for example by modifying its relax- 
ation rate through the Purcell effect [15] . 

In this work, we study an additional consequence of 
the resonator by investigating the driven dynamics and 
the dephasing of a flux qubit [16] that is tunably coupled 
to a harmonic oscillator [T7H20] . We find that the res- 
onator mediates an indirect driving field that interferes 
with the direct drive set by the qubit-antenna coupling, 
thereby modifying both the amplitude and the phase of 
the net driving field. The tunable coupling allows the 
indirect driving to be switched off, but it also opens an 
additional channel for noise to couple into the system. 
Fluctuations in the coupling parameter translate into ef- 
fective driving-field amplitude noise, which reduces the 
qubit coherence during driven evolution. We show that 
the qubit dephasing due to amplitude noise (whether due 
to tunable coupling or otherwise) can be mitigated by a 
rotary echo [21] . a pulse sequence originally developed 
for nuclear magnetic resonance. 

The device, shown in Fig. [TJa), consists of a flux qubit 
and a SQUID embedded in a two- mode LC resonant cir- 



cuit. The diabatic states of the qubit correspond to clock- 
wise or counterclockwise persistent currents in the qubit 
loop [blue arrow in Fig.JlJa)], with energies controlled by 
the flux in the loop. The resonator mode of interest is 
the SQUID plasma mode, depicted by the two red arrows 
in Fig. TJa). The SQUID serves dual purposes: it acts as 
a tunable coupler between the qubit and the resonator, 
and it is also used as a sensitive magnetometer for qubit 
read-out [22]. 

We have investigated two devices with similar lay- 
outs but slightly different parameters, both made of alu- 
minum. Device A was designed and fabricated at MIT 
Lincoln Laboratory and device B was designed and fabri- 
cated at NEC. Figure [TJb) shows a spectroscopy measure- 
ment of device A versus applied flux, with the qubit flux 
detuning <£ qb defined as ^> qb = $ + <I>o/2 and <£ = h/2e. 
The qubit frequency follows / q b = VA 2 + £ 2 , where the 
tunnel coupling A = 2.6 GHz is fixed by fabrication and 
the energy detuning e = 2ip<£> q b//i is controlled by the 
applied flux (Ip is the persistent current in the qubit 
loop). The resonator frequency f r is around 2.3 GHz and 
depends only weakly on <!> q b and 1^. In addition, there 
are features visible at frequencies corresponding to the 
sum of the qubit and resonator frequencies, illustrating 
the coherent coupling between the two systems [3j [23] . 

The system is described by the Jaynes-Cummings 
Hamiltonian [Tl l24ll25] 

H/h = - |cr z + fr (a* a + + y (a + a 1 ") <j z . 

Here, a and a) are the creation/annihilation operators of 
the resonator field and g\ is the dipole coupling between 
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FIG. 1. (a) Circuit diagram of the qubit and the oscillator. 
The qubit state is encoded in currents circulating clockwise 
or counterclockwise in the qubit loop (blue arrow), while the 
mode of the harmonic oscillator is shown by the red arrows, 
(b) Spectrum for device A, showing the qubit and the har- 
monic oscillator. In addition, the two-photon qubit (/ q b/2) 
and the qubit+resonator (/ q b + f T ) transitions are visible, (c) 
Flux induced in the qubit loop by the dc bias current ib- 
The black lines are parabolic fits, (d) First-order coupling 
between the qubit and the ground state of the harmonic os- 
cillator, showing that the coupling is tunable by adjusting ib- 
The coupling is zero at ib = ^) which is slightly offset from 
Ih = due to fabricated junction asymmetry. The derivative 
de/dlb is calculated from the curves in panel (c). The qubit 
parameters are: Ip — 175 n A for device A and Ip — 180 n A 
for device B. The resonators have quality factors Q « 100. 
The right-hand axis is calculated using f T = 2.2 GHz and 
C e ff — 2C— 14 pF for both samples. 



the qubit and the resonator. In this work, we do not 
consider higher-order coupling parameters [24H26] . 

The coupling g\ is mediated by the SQUID. When 
the two SQUID junctions are symmetric, the current of 
the resonator mode splits equally into the two SQUID 
arms, and therefore no net flux is induced into the qubit 
loop. The qubit is thus effectively decoupled from the res- 
onator. However, in the presence of a magnetic field, ap- 
plying a dc bias current Ih creates an asymmetric phase 
drop over the two SQUID junctions. This causes the res- 
onator current to be slightly larger in one of the arms, 
which will produce a flux in the qubit loop. The cou- 
pling to the resonator can thus be controlled in situ by 
changing ib [27] - 

Figure [TJc) shows the flux induced into the qubit loop 
as a function of the dc bias current ib for the two de- 
vices. Given the similarity in the design, both samples 
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FIG. 2. (a) Rabi frequency of qubit A, measured vs ib at 
^ q b = 0. The driving field seen by the qubit contains two 
components: one is due to direct coupling to the antenna, 
the other is due to the coupling mediated by the resonator, 
(b) Rabi traces for a few of the data points in panel (a). 
The microwaves in the antenna have the same amplitude and 
frequency for all traces, (c) Direct coupling between the an- 
tenna and the qubit, extracted from measurements similar 
to the one shown in panel (a). The coupling depends only 
weakly on frequency, (d) Microwave current in the resonator, 
induced by a fixed microwave amplitude in the antenna. The 
black line is a fit to the square root of a Lorentzian, describ- 
ing the oscillation amplitude of a harmonic oscillator with 
f r = 2.3 GHz and Q = 100. 



display similar behavior, with the bias current gener- 
ating a parabolic shift in ^> q b [28l [29] . Since <£ q b con- 
trols the qubit energy detuning e, the first-order qubit- 
resonator coupling strength is determined by the deriva- 
tive de/dl h = (2Jp/ft)(0$ qb /0Ib). The bare coupling 
coefficient between qubit a nd resonato r is then g\ = 
(ds/dIh)SIo, where 5Iq = ^2ii 2 hf^C e ^ is the rms am- 
plitude of the vacuum fluctuations and C e R = 2C is the 
total capacitance of the resonant circuit [25] . The deriva- 
tive de/dlh and the coupling gi are plotted in Fig. [TJd), 
determined at dc from the measured relation between 
$ q b and ib shown in Fig.[TJc). Note that in both devices 
g\ can be tuned over a range of a few tens of MHz, and 
that the coupling is turned off at Ih = The device 
parameters are given in the figure caption. 

Having determined the coupling coefficients, we turn 
to analyzing how the presence of the resonator influences 
the qubit's driven dynamics. Figure a) shows the ex- 
tracted Rabi frequency /R a bi of qubit A as a function 
of ib, measured at / q b = 2.6 GHz. We find that /R a bi 
changes by a factor of five over the range of the mea- 
surement, which is surprising since both the amplitude 
and the frequency of the microwave current I^ enna in 
the antenna are kept constant. The data points were ob- 
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tained by fitting Rabi oscillations to decaying sinusoids, a 
few examples of Rabi traces for different values of Ih are 
shown in Fig. [2^b). We carefully calibrated the flux at 
each data point, to make sure that the qubit was driven 
on resonance and at the degeneracy point (e = 0). 

To explain the results of Fig. |2|a), we need to con- 
sider the indirect driving of the qubit mediated by the 
harmonic oscillator. As seen in Fig.[TJb), at zero flux de- 
tuning the qubit frequency (/ q b = 2.6 GHz) is relatively 
close to the resonator frequency (f r = 2.3 GHz). We 
therefore expect the microwave drive in the antenna to 
off-resonantly induce a microwave current I™ w in the res- 
onator, which is proportional to the square root of the av- 
erage photon population, yj (n). By setting Ih ^ the 



coupling between the resonator and the qubit is turned 
on, and the resonator current I™ w will start driving the 
qubit. To describe this indirect driving, we treat the 
resonator classically and write the qubit Hamiltonian in 
Eq. as 



H qh /h= -1/2 [A<r x + [s' 



dc 



- e mw cos(2tt ft)] a z 



(2) 



(a) 8 



N 

o 

£= 

CD . 




y2 $ 



(b) 



100 



50 



CD 





U 4ntenna ^.2 JlA . 
* (antenna - 8-1 

— Fit, Eqs. (2, 3) 









-5 5 

Flux detuning O qb [mO ] 



-500 



/ b [nA] 



500 



FIG. 3. (a) Spectrum of device B. The spectral line at 
2 GHz is the resonator, whereas the qubit tunnel coupling is 
A = 5.4 GHz. (b) Rabi frequency vs bias current ib, measured 
at / = 5.4 GHz and $ q b = and for two different microwave 
drive amplitudes I^tenna- Similar to device A, the Rabi fre- 
quency depends strongly on lb , and scales lin early with drive 
amplitude. The black lines are fits to Eqs. (3|4), using the 
same coupling parameters for both sets of data. Note that 
the range of lb in Fig. |3jb) is several times larger than in 
Fig- da). 



Here, the drive amplitude e mw experienced by the qubit 
becomes a combination of the drive e 



direct ' 



due to di- 
rect coupling between antenna and qubit, and the drive 
(de/dl\>)l™ w mediated by the resonator. We get: 
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where = 6 a — r is the phase difference between the 
direct drive and the drive mediated by the resonator. 
The Rabi frequency due to the drive e mw depends on the 
qubit's quantization axis, which changes with the static 
energy detuning £dc : 



/Rabi 



2 yii + A^' 



(4) 



Fitting the data in Fig. [2^a) to Eqs. ( |3|4[ ) allows us to 
extract the parameters £^ ct , I™ w and 0. The different 
drive contributions are plotted together with the data in 
Fig.[2|a). The direct drive is independent of ib, while the 
drive I^ w (de/dlh) mediated by the resonator increases 
linearly with \Ih\, which originates from the linear depen- 
dence of gi shown in Fig. [ljd) . The minimum in Rabi 
frequency occurs at a value of lb slightly shifted from 
the point I£ where gi = 0. This offset appears because 
of the phase difference 6 between the two drive compo- 
nents. The fit gives = —75°, which is consistent with a 
resonator driven above its resonance frequency. 

Figures |5|c) and [2^d) show how the two drive com- 
ponents depend on microwave frequency, measured by 
changing the static flux detuning <£ q b to increase the 
qubit frequency [see Fig.[TJb)]. The direct drive only de- 
pends weakly on frequency (due to cable losses), whereas 
the drive mediated by the resonator drops sharply as the 
qubit-resonator detuning increases. The black curve in 



Fig. |2|d) is the frequency response of a harmonic oscil- 
lator with f r = 2.3 GHz and Q = 100, with amplitude 
normalized to match the data. 

To further investigate how the presence of the res- 
onator affects the qubit dynamics at large detunings, we 
performed measurements on device B. Figure [3ja) shows 
a spectrum of that device, where the qubit and the res- 
onator mode (f r = 2 GHz) are clearly visible. This device 
has a larger tunnel coupling (A = 5.4 GHz), which allows 
us to operate the qubit at large frequency detuning from 
the resonator while still staying at e^c — 0, where the 
qubit, to first order, is insensitive to flux noise [28} [29]. 
The qubit-resonator detuning corresponds to several hun- 
dred linewidths of the resonator, which is the regime of 
most interest for quantum information processing [T2] . 

Figure [3](b) shows the Rabi frequency vs bias cur- 
rent Ih of device B, measured at / = 5.4 GHz and 
for two different values of the microwave drive current 
^antenna- Similarly to Fig. [2|a), the Rabi frequency 
clearly changes with /b, but the dependence is weaker 
than in Fig. [2|a) because of the larger frequency de- 
tuning. Note that /R a bi scales linearly with / ^Jenn a 
for all values of Ih,. By fitting the data to Eqs. (3[4), 

We find ^ectffienna = 6.4 MHz/ /iA, CV&na = 

2.4nA//iA and = -155°. The large phase differ- 
ence for device B causes the minimum in /R a bi to shift 
away from the point close to Ih = where the coupling 
#i=0 [see Fig. [ljd)] . We attribute the large phase 
shift to influences from a second resonant mode, which 
is formed by the two L and the two C in the outer loop 
of Fig. [TJa) [I8l [30]. For sample B, this mode resonates 
around 5 GHz. 

The results of Figs. [2] and [3] show that the microwave 
signal mediated by the resonator plays a significant role 
when driving the qubit, appearing already at moderate 
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qubit-resonator coupling gi and persisting even when the 
two systems are far detuned. The design investigated 
here allows the coupling to be turned off [g\ = in 
Fig. [TJd)], but it comes with a drawback: the param- 
eter used to control the coupling (1^ in our setup) also 
provides a way for low-frequency noise to enter the sys- 
tem. Consider the relation between the /R a bi an d lb in 
Fig. |3|b) : fluctuations 51^ near 1^ = will cause fluctua- 
tions in the amplitude of the drive field seen by the qubit, 
which will lead to dephasing during driven evolution. 

To quantify the dephasing, we linearize the relation 
between the Rabi frequency and 1^ close to 1^ = as 
/ /o[l + rSI b ], where f / Ra bi(4 0) and r 
(a/Rabi/5/b)//o = -1.28 (/iA)" 1 is given by Eqs §4) 
or from Fig.[3|b). We model the fluctuations Sly, as nor- 
mally distributed, with standard deviation 07. Assuming 
the noise to be quasi- static, where the value of 5I\> is con- 
stant during a single trial but differs from run to run [31] , 
we find that the Rabi oscillations decay as 
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In addition, the qubit energy relaxation time T\ gives 
an exponential contribution to the Rabi decay, with time 
constant 4Ti/3 given by the Bloch equations. The Rabi 
decay also depends on the flux noise at the Rabi fre- 
quency, but this contribution can be disregarded when 
operating the qubit at £& c = 0, where the qubit is in- 
sensitive to first-order flux noise [32 j. The total decay 
envelope f(t) of the Rabi oscillations becomes 

f(t) = e"*V^ T ^ 2 , with TV = l/(v^7r/ r<7/). (6) 

Note that the Gaussian decay constant due to the 
effective amplitude fluctuations is inversely proportional 
to /o, the average Rabi frequency. This is a consequence 
of having noise in the coupling between the qubit and the 
antenna; the effective amplitude fluctuations seen by the 
qubit will scale with the drive amplitude. 

The red circles in Fig. [4ja) show the envelope of Rabi 
oscillations measured for /R a bi = 65 MHz, together with 
a fit to Eq. The qubit energy relaxation T\ = 11.7 /is 
is known from separate experiments [32], leaving = 
4.3 /is as a fitting parameter. In Fig. Qb), we plot the 
Rabi decay time versus /R a bi? extracted from envelopes 
similar to Fig.[4^a). To capture both the exponential and 
the Gaussian decay, we plot the time T e for the envelope 
to decrease by a factor 1/e. For the lowest Rabi fre- 
quency, the decay time is within 25% of the upper limit 
set by qubit relaxation, but it decreases with /R a bi, as 
expected from Eq. The black solid line shows a fit 
to Eq. giving a value of 07 = 0.8 nA for the noise in 
ib- The effective fluctuations in the drive amplitude are 
rai = 0.06%. We cannot rule out that part of that noise 
may be caused by instrument imperfections. 

Dephasing due to low-frequency fluctuations of the 
qubit frequency is routinely reduced by performing a 
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FIG. 4. (a) Decay envelopes of the Rabi and rotary-echo 
sequences for device B, measured with /R a bi = 65 MHz. The 
solid lines are fits to Eq. (b) Decay times for Rabi and 
rotary echo, extracted from fits similar to the ones shown in 
panel (a). The dashed line shows the upper limit set by qubit 
energy relaxation. The dotted line marks the position for the 
decay envelope shown in panel (a). (c,d) Schematic diagrams 
describing the two pulse sequences in (a) and (b). For rotary 
echo, the phase of the microwaves is rotated by 180° during 
the second half of the sequence. 



Hahn-echo experiment [33]. Similarly, the fluctuations 
in drive field that cause decay of the Rabi oscillations 
in Fig. Qb) can be mitigated with the rotary-echo pulse 
sequence [21], depicted in Fig.[4jc), which for driven sys- 
tems is analogous to the Hahn-echo sequence. By shifting 
the phase of the drive by 180° after a time t p /2, any addi- 
tional rotations, acquired due to slow fluctuations in the 
drive amplitude during the first half of the sequence, will 
cancel out during the reversed rotations in the second 
half of the sequence. 

The blue squares in Fig. [4ja) shows the decay of the 
rotary-echo sequence, measured for /Rabi = 65 MHz. The 
rotary-echo data shows a clear improvement compared 
to the Rabi decay for the same parameters [red circles in 
Fig.Qa)]. We fit the rotary-echo data to Eq. (|6| and plot 
the extracted decay times together with the results from 
Rabi measurements in Fig. Qb). The rotary-echo signal 
outperforms the Rabi decay over the full range of Rabi 
frequencies, and reaches the upper limit set by qubit re- 
laxation (4Ti/3) at low frequencies. For intermediate fre- 
quencies, the rotary-echo decay times are slightly shorter 
than 4Ti/3; we attribute the reduced coherence times to 
fluctuations in 1^ that occur on time scales comparable 
to the length of the pulse sequence. Noise at frequen- 
cies around l/t v will not be refocused by the reversed 
drive pulse, since the rotary-echo sequence has similar 
filtering properties as the Hahn-echo [21] . At the highest 
drive amplitudes (/R a bi > 100 MHz), we observe a strong 
increase in decoherence, probably due to heating. The 
indirect driving can also be reduced by driving the qubit 
with two antennas with different amplitudes and phases 
[34] , but it requires a more complicated setup. 

To conclude, we have investigated interference effects 
occurring when driving a qubit that is tunably coupled 
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to a harmonic oscillator. Although the addition of a 
coupling control parameter opens up an extra channel 
for dephasing, we show that its influence is reversible 
with dynamical decoupling techniques. The results are 
relevant for any type of qubit that is tunably coupled 
to a resonator, and they show that despite engineering 
limitations, imperfections can be reversed by applying 



proper decoupling protocols. In analogy with multi-pulse 
Hahn-echo experiments, we expect the incorporation of 
additional rotary echos to further improve the coherence 
times [32] . 
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